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<n : 

£h ! Abstract 

a . 

We discuss the possibility that in finite density QCD an anisotropic 
phase is realized. This case might arise for quarks with different chem- 
ical potential and/or different masses. In this phase crystalline struc- 
t— I ' tures may be formed. We consider this possibility and we describe, in 

the context of an effective lagrangian, the corresponding phonons as 
\^ | the Nambu-Goldstone bosons associated to the breaking of the space 

I ' symmetries. 

o : 

1 Introduction 

Ph! The study of color superconductivity goes back to the first days of QCD [1], 

o ,| but only recently this phenomenon has received a lot of attention (for recent 

reviews see refs. [2, 3]). Naively one would expects that, due to asymptotic 
freedom, quarks at very high density would form a Fermi sphere of almost free 
fermions. However, Bardeen, Cooper and Schrieffer [4] proved that the Fermi 
surface of free fermions is unstable in presence of an attractive, arbitrary 
small, interaction. Since in QCD the gluon exchange in the 3 channel is 
attractive one expects the formation of a coherent state of particle/hole pairs 
(Cooper pairs). An easy way to understand the origin of this instability is 
to remember that, for free fermions, the Fermi energy distribution at zero 
temperature is given by f{E) = 6(fi — E) and therefore the maximum value 
of the energy (Fermi energy) is Ep = fi. Then consider the corresponding 
grand-potential, Q = E — /iiV, with /i = Ep. Adding or subtracting a 
particle (or adding a hole) to the Fermi surface does not change Q, since 
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Q — > (E ± — //(TV =b 1) = Q. We see that the Fermi sphere of free 
fermions is highly degenerate. This is the origin of the instability, because if 
we compare the grand-potential for adding two free particles or two particles 
bounded with a binding energy Eb, we find that the difference is given by 
il B — Q = —E B < 0. Since a bound state at the Fermi surface can be 
formed by an arbitrary small attractive interaction [5], it is energetically 
more favorable for fermions to pair and form condensates. 

From the previous considerations it is easy to understand why to realize 
superconductivity in ordinary matter is a difficult job. In fact, one needs an 
attractive interaction to overcome the repulsive Coulomb interaction among 
electrons, as for instance the one originating from phonon exchange for elec- 
trons in metals. On the other hand the interaction among two quarks in the 
channel 3 is attractive , making color superconductivity a very robust phe- 
nomenon. Notice also that once taken into account the condensation effects, 
at very high density one can use asymptotic freedom to get exact results. 
For instance, it is possible to get an analytical expression for the gap [6]. In 
the asymptotic regime it is also possible to understand the structure of the 
condensates. In fact, consider the matrix element 

(oiC^lo) (i) 

where a, (3 = 1,2,3 are color indices, a,b — 1,2 are spin indices and i,j = 
1, • • • ,N are flavor indices. Its color, spin and flavor structure is completely 
fixed by the following considerations: 

• antisymmetry in color indices (cx,/3) in order to have attraction; 

• antisymmetry in spin indices (a, b) in order to get a spin zero conden- 
sate. The isotropic structure of the condensate is favored since it allows 
a better use of the Fermi surface; 

• given the structure in color and spin, Pauli principles requires antisym- 
metry in flavor indices. 

Since the momenta in a Cooper pair are opposite, as the spins of the quarks 
(the condensate has spin 0), it follows that the left(right)-handed quarks can 
pair only with left(right)-handed quarks. In the case of 3 flavors the favored 
condensate is 

3 

(o|Cl^|o> = -(o\€r^r\°) = A E eQ/3C W (2) 

c=\ 
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This gives rise to the so-called color-flavor-locked (CFL) phase [7, 8]. The 
reason for the name is that simultaneous transformations in color and in 
flavor leave the condensate invariant. In fact, the symmetry breaking pattern 
turns out to be 

SU(3) C <g) SU(3) L ® SU(3) R <g> U(l) B -> SU(3) C+L+R 

where SU(3) c +l+r is the diagonal subgroup of the three SU(3) groups. This 
is the typical situation when the chemical potential is much bigger than the 
quark masses m u , raj and m s (here the masses to be considered are in prin- 
ciple density depending). However we may ask what happens decreasing the 
chemical potential. At intermediate densities we have no more the support 
of asymptotic freedom, but all the model calculations show that one still has 
a sizeable color condensation. In particular if the chemical potential /i is 
much less than the strange quark mass one expects that the strange quark 
decouples, and the corresponding condensate should be 

(OICl^IO) = Ae*% (3) 

In fact, due to the antisymmetry in color the condensate must necessarily 
choose a direction in color space. Notice that now the symmetry breaking 
pattern is completely different from the three-flavor case. In fact, we have 

SU(3) C ® SU(2) L ® SU(2) R ® U{1) b -> SU(2) C ® SU(2) L ® SU(2) R ® U(l) B 

It is natural to ask what happens in the intermediate region of /i. It turns 
out that the interesting case is for /i M s 2 /A. To understand this point let 
us consider the case of two fermions: one massive, mi = M s and the other 
one massless, m 2 = 0, at the same chemical potential \x. The Fermi momenta 
are of course different 

p Fl = vV - M?, Pf 2 = A* (4) 

The grand potential for the two unpaired fermions is (factor 2 from the spin 
degrees of freedom) 

SW. = 2 p ^gj (v^T^ - „) + 2 j'" dp I - ,) (5) 
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In order to pair the two fermions must reach some common momentum p^ mm , 
and the corresponding grand potential is 

ffcomm rJ^Tl f I \ P Pcomm /7^T) 



o 

2 A2 



4tt 2 



(6) 



where the last term is the energy necessary for the condensation of a fermion 
pair [9]. The common momentum Pc 0mm can be determined by minimizing 
^pair. with respect to Pc 0mm , with the result 

f 

Pcomm = V ~ (7) 

It is now easy to evaluate the difference n unpair . — fi pa i r . at the order Mf, with 
the result 

^unpair. - ^pair. « (M* ~ 4AV ) (8) 

We see that in order to have condensation the condition 

M 2 

must be realized. The problem of one massless and one massive flavor has 
been studied in ref. [10]. However, one can simulate this situation by taking 
two massless quarks with different chemical potentials, which is equivalent 
to have two different Fermi spheres. The big advantage here is that one can 
use a study made by Larkin and Ovchinnikov [11] and Fulde and Ferrel [12]. 
These authors studied the case of a ferromagnetic alloy with paramagnetic 
impurities. The impurities produce a magnetic field which, acting upon the 
electron spins, gives rise to a different chemical potential for the two popu- 
lations of electrons. It turns out that it might be energetically favorable to 
pair fermions which are close to their respective Fermi surface (LOFF phase). 
However, since the Fermi momenta are different, the Cooper pair cannot have 
zero momentum and there is a breaking of translational and rotational invari- 
ance. Therefore, a crystalline phase can be formed. The previous situation 
is very difficult to be realized experimentally, but there have been claims 
of observation of this phase in heavy-fermion superconductors [13] and in 
quasi-two dimensional layered organic superconductors [14]. The authors of 
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ref. [15] have extended the calculation of ref. [11, 12] to the case of two-flavor 
QCD and we will review here their results. Also, since the LOFF phase can 
give rise to crystalline structures, phonons are expected. We will also discuss 
the effective lagrangians for the phonons in different crystalline phases and 
show how to evaluate the parameters characterizing them, as the velocity of 
propagation. Finally we will consider an astrophysical application. 



According to the authors of ref. [11, 12] when fermions belong to two different 
Fermi spheres, they may prefer to pair staying as much as possible close to 
their own Fermi surface. When they are sitting exactly at the surface, the 
pairing is as shown in Fig. 1. We see that the total momentum of the pair 



Figure 1: Pairing of fermions belonging to two Fermi spheres of different 
radii according to LOFF. 

is pi + f>2 = 2<f and, as we shall see, \q \ is fixed variationally whereas the 
direction of q is chosen spontaneously. Since the total momentum of the pair 
is not zero the condensate breaks rotational and translational invariance. 
The simplest form of the condensate compatible with this breaking is just a 
simple plane wave (more complicated functions will be considered later) 



It should also be noticed that the pairs use much less of the Fermi surface than 
they do in the BCS case. In fact, in the case considered in Fig. 1 the fermions 
can pair only if they belong to the circles depicted there. More generally there 



2 The LOFF phase 




(10) 
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is a quite large region in momentum space (the so called blocking region) 
which is excluded from the pairing. This leads to a condensate smaller than 
the BCS one. 

Before discussing the LOFF case let us review the gap equation for the 
BCS condensate. We have said that the condensation phenomenon is the 
key feature of a degenerate Fermi gas with attractive interactions. Once one 
takes into account the condensation the physics can be described using the 
Landau's idea of quasi-particles. In this context quasi-particles are nothing 
but fermionic excitations around the Fermi surface described by the following 
dispersion relation 



e(p, A BCS ) = y/? + A% cs (11) 

with 

Z = E(p)-im?^ -(p-f F )=H F -(f-p F ) (12) 

Op p=p F 

The quantities v F and {p — p F ) are called the Fermi velocity and the resid- 
ual momentum respectively. A easy way to understand how the concept of 
quasi-particles comes about in this context is to study the gap equation at 
finite temperature. For simplicity let us consider the case of a four-fermi 
interaction. The euclidean gap equation is given by 

19 J (27r)MP4-^) 2 +|p| 2 + A| C5 1 } 

From this expression is easy to get the gap equation at finite temperature. 
We need only to convert the integral over p 4 into a sum over the Matsubara 
frequencies 

r d 3 p 1 

1 = gT J ^ n Ep+W+^P^) (14) 

Performing the sum we get 

g f d 3 p 1 - n u - n d 



II 



(2tt)3 e(p,A BCS ) 



(15) 



Here n u and are the finite-temperature distribution functions for the ex- 
citations (quasi-particles) corresponding to the original pairing fermions 

n u = ^= e ^ AB a s)/ T + I ( 16 ) 
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At zero temperature (n u = n d — > 0) we find (restricting the integration to a 
shell around the Fermi surface) 

£ /" dQ p p 2 F d£ 1 

In the limit of weak coupling we get 

A S c5«2£"e- 2 /^ (18) 

where £ is a cutoff and 

P=#- (19) 

IT Dp 

is the density of states at the Fermi surface. This shows that decreasing the 
density of the states the condensate decreases exponentially. 

Let us now consider the LOFF case. For two fermions at different densities 
we have an extra term in the hamiltonian which can be written as 

Hj = -5fia 3 (20) 

where, in the original LOFF papers [11, 12] 5/j, is proportional to the magnetic 
field due to the impurities, whereas in the actual case 5[i = (/ii — fJ>2)/2 and 
(T3 is a Pauli matrix acting on the space of the two fermions. According to 
ref. [11, 12] this favors the formation of pairs with momenta 

Pi = k + q, p 2 = -k + q (21) 

We will discuss in detail the case of a single plane wave (see eq. (10)). The 
interaction term of eq. (20) gives rise to a shift in £ (see eq. (12)) due both 
to the non-zero momentum of the pair and to the different chemical potential 

S = E{0)-n^ E{±k + $-n^5n*i^ii (22) 

with 

Jx = 5/i — vp ■ q (23) 

Here we have assumed 5/i \i (with \i = (/xi + ^2)/ 2) allowing us to expand 
E at the first order in q (see Fig. 1). The gap equation has the same formal 
expression as in eq. (15) for the BCS case 



d 3 p l-n u -n d ^ 



(2vr)3 e(p,A) 
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but now n u ^ n d 

Uu > d = e {e(pA)±fi)/T + i ( 25 ) 
where A is the LOFF gap. In the limit of zero temperature we obtain 

g f d 3 p 



2 



/(^^A)( 1 -"(- e -^-"(- £ + ^ < 26 > 



The two step functions can be interpreted saying that at zero temperature 
there is no pairing when e(p, A) < This inequality defines the so called 
blocking region. The effect is to inhibit part of the Fermi surface to the 
pairing giving rise a to a smaller condensate with respect to the BCS case 
where all the surface is used. 

We are now in the position to show that increasing 5 fx from zero we have 
first the BCS phase. Then at 5 fx = 5fx\ there is a first order transition to the 
LOFF phase [11, 15], and at Sfx = 5fx 2 > there is a second order phase 
transition to the normal phase (with zero gap) [11, 15]. We start comparing 
the grand potential in the BCS phase to the one in the normal phase. Their 
difference is given by 

2 

^BCS — ^normal = ~ i^BCS ~ 25/i 2 ) (27) 

where the first term comes from the energy necessary to the BCS conden- 
sation (compare with eq. (6)), whereas the last term arises from the grand 
potential of two free fermions with different chemical potential. We recall 
also that for massless fermions p F = fx and v F — 1. We have again assumed 
5/x fi. This implies that there should be a first order phase transition from 
the BCS to the normal phase at 5ji = Abcs / V%, since the BCS gap does 
not depend on Six. The situation is depicted in Fig. 2. In order to compare 
with the LOFF phase we will now expand the gap equation around the point 
A = (Ginzburg-Landau expansion) in order to explore the possibility of a 
second order phase transition. Using the gap equation for the BCS phase in 
the first term on the right-hand side of eq. (26) and integrating the other 
two terms in £ we get 

«* log ^ = jvLf^ inh m (28) 



2tt 2 v f & A 2n 2 v F J 4vr A 

where 

C{&) = a/ (5li - qv F cos 6) 2 - A 2 (29) 
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Figure 2: The grand potential and the condensate of the BCS and LOFF 
phases vs. SfM. 

For A — ► we get easily 

,og^ = i/ +1 lo g fl-^, z = lE (30) 



28/1 2 J_ ± V z ) qvp 

This expression is valid for 8\x smaller than the value 5 1x2 at which A = 0, 
therefore the right-hand side must reach a minimum at 8fi = 8fi2- The 
minimum is fixed by the condition 

-tanh- = l (31) 

z z 

implying 

qv F fn 1.2 5n (32) 
Putting this value back in eq. (30) we obtain 

(5/i 2 ^ 0.754 A BCS (33) 

From the expansion of the gap equation around 8^2 it is easy to obtain 

A 2 « 1.76 8fi 2 (8fi 2 -5(i) (34) 

According to ref. [9] the difference between the grand potential in the super- 
conducting state and in the normal state is given by 

^LOFF — ^normal = ~ [ ^f^ 2 (35) 

JO 9 
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Using eq. (18) and eq.(33) we can write 



dg_ = p (IAbcs = p d5fi 2 , . 

~ 2 A BCS ~ 2 8 fa [ ' 

Noticing that A is zero for 5p 2 = 8p we are now able to perform the integral 
(35) obtaining 

^LOFF - ^normal ~ -0.44 p(5fl - 5 faf (37) 

We see that in the window between the intersection of the BCS curve and 
the LOFF curve in Fig. 2 and 8 fa the LOFF phase is favored. Furthermore 
at the intersection there is a first order transition between the LOFF and 
the BCS phase. Notice that since 8 fa is very close to 8 fa the intersection 
point is practically given by 8 fa. In Fig. 2 we show also the behaviour of the 
condensates. Altough the window (8 fa, 8 fa) ~ (0.707, 0.754) A BC s is rather 
narrow, there are indications that considering the realistic case of QCD [16] 
the windows may open up. Also, for different structures than the single plane 
wave there is the possibility that the windows opens up [16]. 



3 Crystalline structures 

The ground state in the LOFF phase is a superposition of states with different 
occupation numbers (N even) 

\0)loff = 

E c ^> ( 38 ) 

N 

Therefore the general structure of the condensate in the LOFF ground state 
will be 

mx)iP(x)) =Y,CNCN + 2e W (N\ii{x)i,{x)\N + 2) = ^A,e*' f (39) 

N N 

The case considered previously corresponds to all the Cooper pairs having the 
same total momentum 2q. A more general situation, although not the most 
general, is when the vectors reduce to a set defining a regular crystalline 
structure. The corresponding coefficients (linear combinations of subsets 
of the Aat's) do not depend on the vectors <^ since all the vectors belong to 



10 



the same orbit of the group. Furthermore all the vectors have the same 
lenght [17] given by eq. (32). In this case 

(0\i/>(x)i/>(x)\Q) = A^e**'* (40) 

i 

This more general case has been considered in [11, 17] by evaluating the 
grand-potential of various crystalline structures through a Ginzburg-Landau 
expansion, up to sixth order in the gap [17] 

Q = aA 2 + ^A 4 + |A 6 (41) 

These coefficients can be evaluated microscopically for each given crystalline 
structure. The procedure that the authors of ref. [17] have followed is to 
start from the gap equation represented graphically in Fig. 3. Then, they 

^ o 



Figure 3: Gap equation in graphical form. The thick line is the exact prop- 
agator. The black dot the gap insertion. 

expand the exact propagator in a series of the gap insertions as given in 
Fig. 4. Inserting this expression back into the gap equation one gets the 

Figure 4: The expansion of the propagator in graphical form. The darker 
boxes represent a A* insertion, the lighter ones a A insertion. 

expansion illustrated in Fig. 5. On the other hand the gap equation is 
obtained minimizing the grand-potential (41), i.e. 

aA + (3 A 3 + 7 A 5 + • • • = (42) 

Comparing this expression with the result of Fig. 5 one is able to derive the 
coefficients a, [3 and 7. 
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Figure 5: The expansion of the gap equation in graphical form. Notations 
as in Fig. 4. 

In ref. [17] more than 20 crystalline structures have been considered, 
evaluating for each of them the coefficients of Eq. (41). The result of this 
analysis is that the face-centered cube appears to be the favored structure 
among the ones considered (for more details see ref. [17]). 



4 Phonons 

Since in the LOFF phase translational and rotational symmetries are broken, 
we expect the corresponding Nambu-Goldstone bosons (phonons) to appear 
in the theory. The number and the features of the phonons depend on the 
particular crystalline structure. We will consider here the case of the single 
plane- wave [18] and of the face-centered cube [19]. We will introduce the 
phonons as it is usual for NG bosons [18], that is as the phases of the con- 
densate. Considering the case of a single plane-wave we introduce a scalar 
field Q(x) through the replacement 

A(x) = exp 2 *^ A -> e^A (43) 

We require that the scalar field <&(x) acquires the following expectation value 
in the ground state 

($(x)) = 2q-x (44) 

The phonon field is defined as 

—<j)(x) = $(x) -2q-x (45) 

Notice that the phonon field transforms nontrivially under rotations and 
translations. From this it follows that non derivative terms for 4>(x) are not 
allowed. One starts with the most general invariant lagrangian for the field 
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<&(#) in the low-energy limit. This cuts the expansion of $ to the second 
order in the time derivative. However one may have an arbitrary number of 
space derivative, since from eq. (44) it follows that the space derivatives do 
not need to be small. Therefore 

£ P ho„o„ = y ^ 2 + E c ^(V 2 ) fc $ j (46) 

Using the definition (45) and keeping the space derivative up to the second 
order (we can make this assumption for the phonon field) we find 

Vonon = g " • V ± - ^V||0 • V||0) (47) 

where 

— * 

V|| = ra(n- V), V ± = V-V||, n = jL (48) 

We see that the propagation of the phonon in the crystalline medium is 
anisotropic. 

The same kind of considerations can be made in the case of the cube. The 
cube is defined by 6 vectors <fj pointing from the origin of the coordinates to 
the vertices of the cube parameterized as in Fig. 6. 

The condensate is given by [17] 

8 3 

Af^A^e^^A (? me%x% (49) 

k=l i=l,(ei=±) 

We introduce now three scalar fields such that 

= 2\q \xi (50) 

through the substitution 

3 

A{x)^A e iei * ( °W (51) 

i=l,(ei=±) 



and the phonon fields 



\<l)^(x) = ^(x)-2\q\x i (52) 
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(-1,-1,1) (1,-1,1) 




(-1,1,-1) (1,1,-1) 



Figure 6: The figure shows the vertices and corresponding coordinates of the 
cube described in the text. Also shown are the symmetry axis. 

Notice that the expression (51) is invariant under the symmetry group of the 
cube acting upon the scalar fields This group has three invariants 

for the vector representation 

I 2 {X) = \X | 2 , h{X) = X\X\ + XlX\ + XjXl h{X) = X\X\Xl (53) 

Therefore the most general invariant lagrangian invariant under rotations, 
translations and the symmetry group of the cube, at the lowest order in the 
time derivative, is 

Vonon = y E (^ W ) 2 + ^s(/2(V$«),/4(V<l>»),/ 6 (V<l>«)) (54) 
1=1,2,3 

Expanding this expression at the lowest order in the space derivatives of the 
phonon fields one finds [19] 

w = \ E \ E \^^\ 2 - \ E w } ) 2 

i=l,2,3 8=1,2,3 i=l,2,3 

- c E di^djt® (55) 

i<i=l,2,3 

The parameters appearing in the phonon lagrangian can be evaluated 
following the strategy outlined in [20, 21]. One starts from the QCD la- 
grangian and derives an effective lagrangian describing fermions close to the 
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Fermi surface, that is at momenta such that p ~ pp but p — pp ^> A. The 
relevant degrees of freedom are the fermions dressed by the interaction, the 
so-called quasi-particles [22]. Going closer to the Fermi surface the gapped 
quasi-particles decouple and one is left with the light modes as NG bosons, 
phonons and un-gapped fermions. It is possible to derive the parameters 
of the last description by the one in terms of quasi-particles evaluating the 
self-energy of the phonons (or the NG bosons) through one-loop diagrams 
due to fermion pairs. The couplings of the phonons to the fermions are ob- 
tained noticing that the gap acts as a Majorana mass for the quasi-particles. 
Therefore the couplings originate from the substitutions (43) and (51). In 
this way one finds the following results: for the single plane-wave 



5 Astrophysical consequences 

A typical phenomenon of the pulsars are the glitches (for a review see [23]), 
that is sudden jumps in the period of the star. If pulsars are neutron stars 
with a dense metallic crust, the effect is explained assuming that some an- 
gular momentum is stored in the vortices present in the inner neutron su- 
perffuid. When the period of the star slows down due to the gravitational 
radiation, the vortices, which are pinned to the crystalline crust, do not par- 
ticipate in the slowing-down until they become unstable releasing suddenly 
the angular momentum. Since the density in the inner of a star is a function 
of the radius, it results that one has a sort of laboratory to study the phase 
diagram of QCD at zero temperature, at least in the corresponding range 
of densities. A possibility is that one has a CFL state as a core of the star, 
then a shell in the LOFF state and eventually the exterior part made up of 
neutrons. Since in the CFL state the baryionic number is broken there is su- 
perfluidity. Therefore the same mechanism explained above might work with 
vortices in the CFL state pinned to the LOFF crystal. This could reinforce 
the ideas about the existence of strange stars (made of up, down and strange 




(56) 



and for the cube 



a 




(57) 



quarks) . 
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